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1 Introduction 

The entanglement entropy is a useful universal quantity when we would like to employ the AdS/CFT 

correspondence [1] to study non-perturbative aspects of string theory as quantum gravity. In [2], 

it has been proposed that the entanglement entropy in conformal field theories (CFTs) can be 

calculated from the area of minimal surface in anti de-Sitter (AdS) spaces. It is explicitly given by 

the formula 

Area(7A) , . 

where Gm is the Newton constant in the Einstein gravity on the AdS space; Sa is the entanglement 
entropy for the subsystem A, which can be chosen arbitrarily; ^^ is the codimension two minimal 
area surface which ends on the boundary of A. Also we require ^a is homologous to A. Even 
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though there has been no precise proof of this formula ([T]), many non-trivial evidences have been 
accumulated until now. For example, the strong subadditivity, which is one of the most important 
inequality satisfied by the entanglement entropy, has been shown in [3]. Quite recently, the paper [1] 
showed that all known inequalities are satisfied by the holographic entanglement entropy (HEE) 
given by ([T|) and moreover HEE leads to more constrained inequalities than entanglement entropies 
in general quantum field theories. Also, the agreements of the logarithmic terms in both side has 
been shown in [2} [51-110] . Moreover, in [11], a highly non-trivial consistency check has been made 
when the subsystem A is two intervals [12] in two dimensional CFTs by employing so called Renyi 
entropies. Refer to [13] and references therein for more evidences. For quantum field theoretic 
analysis of entanglement entropy refer e.g. to the excellent reviews [14]I15]. 

The holographic entanglement entropy ([T|) assumes that the gravity on the AdS space is defined 
by the Einstein gravity with a negative cosmological constant. In the CFT side, this means that 
we are taking the large 't Hooft coupling limit. To analyze a full quantum gravity which appears 
in string theory and to extend the results to the weak coupling region of the CFT dual, we need 
to include quantum corrections which are effectively described by higher derivative terms. In the 
recent papers [16[I17] studied the corrections of holographic entanglement entropy in the presence 
of higher derivative terms (see also [18] for related discussions). Though it is very complicated 
to find a holographic entanglement entropy formula for general higher derivative corrections as 
noted in |16] . a class of theories called Lovelock gravities, the authors of [161117] found important 
evidences that natural extensions of ([T]) give the correct holographic entanglement entropy. In the 
Ti? Gauss-Bonnet gravity, which is the simplest example of Lovelock gravities, it is given by the 
following formula (this expression itself has been already speculated in [19] ) 

^^=4^/ ^(1 + ^^'^)' (2) 

where rj is the coupling constant of Gauss-Bonnet term, which will be explained later; h is the 
induced metric on the three dimensional surface ja chosen arbitrarily; TZ is the intrinsic curvature 
of h. We again require that 'ja has a boundary which coincides with that of A and we need to add 
a boundary Gibbons-Hawking term to ([2]) in actual calculations. The profile of 7a is determined 
by minimizing the functional ([2]). 

The purpose of this paper is to investigate further the properties of higher derivative corrections 
of the entanglement entropy. In the first half of this paper, we will consider explicit setups of 
string/M-theory and calculate the higher derivative corrections of the holographic entanglement 
entropy. Since the leading higher derivative correction is the quartic polynomial of the Weyl 
curvature [20]I21]. this is not in the class of Lovelock gravities. To obtain non-trivial results in a 
tractable way, we will concentrate on the examples of AdS soliton spaces, which are dual to confining 



gauge theories [55], and choose the subsystem as just the half of the total space. In the type IIB 
string case, this is dual to the subleading corrections to the entanglement entropy in the A/" = 4 
super Yang-Mills (SYM) on a circle with the anti-periodic boundary condition for fermions. In the 
latter half of this paper, we will study the holographic entanglement entropy in the 4D and 5D 
Gauss-Bonnet gravity, especially when the subsystem A is defined by a strip. We will analyze the 
behavior of the entropy when we change the width of the strip and see how the phase transition is 
affected by the higher derivatives. The holographic entanglement entropy for AdS solitons without 
higher derivative corrections has been already studied in |23fl25] . Later, qualitatively similar results 
has been recently reproduced in lattice gauge theory approaches |26H28| . 

This paper is organized as follows: In section [2] we calculate the higher derivative corrections of 
holographic entanglement entropy for the AdSs soliton in IIB string and for the AdS4j soliton in 
M-theory. In section [3l we compute the holographic entanglement entropy for AdS soliton in the 
Gauss-Bonnet gravity in four and five dimension. We choose the subsystem ^ to be a half space 
in all of the above calculations. In section HI we analyze the holographic entanglement entropy for 
AdS soliton in the five dimensional Gauss-Bonnet gravity when A is give by a strip. In section [5l 
we summarize our conclusions. 

2 Holographic Entanglement Entropy with Higher Derivative Cor- 
rections in String and M Theories 

In this section, we will investigate the leading higher derivative corrections to the holographic 
entanglement entropy in string theory and M-theory. We consider AdSs soliton geometry in type 
IIB superstring and AdS4 and AdSy soliton geometries in M-theory. We choose the subsystem A 
as the half of the total space. 

2.1 Type IIB Superstring 

Consider the AdSsxS^ solution in type IIB string theory. The relevant part of the type IIB 
supergravity action with the leading a' correction looks like in the Euclidean signature |20j 



IGvrG^ 



h^l'^""^ 



n-^{d,cpf-^^Fi + ... + ^.e-l'f>w 



(3) 



where we keep the relevant terms which are linear with respect to 7, which is supposed to be very 
small. There are also other such terms depending on other antisymmetric fields, represented by 
dots, but they do not have any contributions to our discussions below. W is defined by the following 



7^^ term by using the Weyl curvature C 






and the constant 7 is give by 



7 = k(3)«"- (5) 



The AdSsX S^ is the solution to this theory with the higher derivatives. This type IIB string 
background is dual to the four dimensional AA = 4 SU{N) super Yang- Mills theory [1]. The 
standard dictionary tells us that the AdS radius L is related to the 't Hooft coupling A = 2-KgsN 

by 

L^ = 4TTgsa''^N = 2X0'"^. (6) 

Therefore taking into account the linear order of 7 is to consider the next leading order correction 
in the strong coupling expansion of A. The Yang-Mills coupling is related to the string coupling 

9ym = '^'^9s- (7) 



The ten dimensional Newton constant is given by 

4^2 

2iV2- 

2.1.1 AdSs soliton and the large 't Hooft coupling limit of TV = 4 SYM 



G'P = 8n^a'-gt = ^. (8) 



Consider the AdSs soliton background, which is dual to the A/" = 4 SYM compactified on a circle S^ 
(radius R) with the anti-periodic boundary condition for fermions [22]. Its leading order solution 
(i.e. 7 = 0) is given by 

ds' = ^[dtl + dx' + dy'] + ^Y^p + ^(1 - r^o/r')de' + L^d^l (9) 

The periodicity of 9 is given by ^ ~ + 2TrR with 

R = ^. (10) 

2ro 

We define the subsystem A by the following half space on the time slice tE = 

x>0, -oo<y<oo, 0<e<2TrR. (11) 

The contribution of holographic entanglement entropy in the leading order of 7 is given by the area 
law formula dTJ. In the AdS soliton case, we obtain explicitly [23] 

(0) _ Area(72^ _ ^3l4 ^L2 
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where we employed 
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(5) 

N 



G 



(10) 

N 



■kL^ 



vol(55) 2iV2' 



(13) 



using vol (5^ 



3r5 



vrU 



2.1.2 Higher derivative corrections to the entanglement entropy 



Next we would like to calculate the higher derivative corrections to the previous calculation (J12p . 
To calculate the holographic entanglement entropy, we place a deficit angle 6 = 27r(l — n) on the 
three dimensional surface 7a defined hy x = tE = [2l[T9l[29] . When the extrinsic curvature of the 
submanifold ja is vanishing, the curvature tensor behaves like (assuming 1 — n is infinitesimally 
small) [30] 






Nf,^N^p)6{^, 



(14) 



where we defined N^y = n^'^n^jf +n^n^. The unit vectors n|f and n^ are orthogonal to the surface 
7A- T^liu etc. denotes the original curvature tensor without the deficit angle. 



By plugging p^ into ([3]), we obtaiqf 

W = ^^^ + 27r(n - 1)5{ja] 



^8^16 ^ ^-v'^^ L6r 

Therefore, the contribution from the ly-term is estimated to be 



4nr-12 



S 



(W) 



d_ 

dn 



log 



7 



7^ 



zrT 



levrG 



(10) 

N 



dx^ 



n=l 



'9lA 



lA 



807rr£"l 

^6^12 






VG 



(5) 



where ^/g^ 



7A^Vi 
rL'^dQ^ is the induced volume factor on 7^4. 



(15) 



(16) 



Next we need to consider the higher derivative corrections of the metric. It is given by 



ds' 



e-^^dsi + e'^idn^r, 



(17) 



By naive substitution, we have several terms including ^(7^)'' [k > 2) in the subleading part 0((1 



Precisely speaking, we should regularize the (5-functions at first and take the limit in the last stage of the calculation. 
But anyway this does not affect our results below. 



where v is order 0(7). The five dimensional part is given by a modification of the AdS soliton 
metric [3n[32l 



where 



ds^ = f(r)de^ + -^dr^ + r^idtl + dx'^ + dy^), (18) 

9{r) 



fir) = '-^^ (1 + 1 -air)), 

gir) = '-^^{l+7-m)- (19) 



The functions a(r) and /3(r) are defined by 



15 / r^ r« r^2\ 



15 / r,^ -» -12 



The radius in the ^-direction can be found by requiring the smoothness of the metric 

2.fi = ^-l^ + 0h'). (21) 

Since we aheady took care of the W term in ()16p , here we only consider the contribution from 
the Einstein-Hilbert term. This is simply given by the area law formula as 

^2 ^2 ^2 ■ 



^{EH) _ Vi 27rR fri^ r^ lo^^M 



'N 



-r^ 



2L4 °° SR ARL^ 

Finally the total expression of Sa up to 0(7) is given by 

2L4 '^~ 8R ARL^ 

r2i 



+ ^^^^7^ + 0(7^). (22) 



^'^1 ^ 1 , C(3),,,,_3/2 



- (d-) + Vl-8 + ^(^^)'^^^J' ^^^) 

where div. denotes the quadratically divergent term. We are interested in its finite term, which is 
independent of the UV cut off. This final result means that the entropy increases as the 't Hooft 
coupling gets decreased. Indeed, the free field theory result obtained in [23] is given by 



-.free . ,. n N Vi 

'^ =(^^"-)-T2¥ 

and matches with this expectation from our gravity calculation. 



ST' = (div.) - ^, (24) 



2.2 M-Theory 

It will be interesting to compute the higher derivative corrections of the holographic entanglement 
entropy in M-theory backgrounds. The one- loop corrected action of the IID supergravity is given 
by [21] 

I = -^jdx''^[n+... + j-W + ...], (25) 

where the dots represent the terms depending on antisymmetric forms again. W is given by the 

same form as (j3]), and 7 is 

47r2 1 
7 = — ^i\- (26) 

The IID Newtonian constant G^ and the Planck length £p are related to kh as 

SvrGii') = KJ^ = 2^7r^i%. (27) 

2.2.1 M2-brane Soliton 

Consider the AdS4 soliton which is dual to a three dimensional CFT on N M2-branes compactified 
on a circle. First of all, upon dimensional reduction on an 5^, the correction term /(^) in (|25p is 
reduced tec 

I^"^^ = ^ fd^xV9~47W,, (28) 



levrGg) 



(4) 



where W4 is given by the same form as Q in 4D and the 4D Newtonian constant Gj^ is given by 

^(11) 

< = ^WY- (29) 

voI(d') 

Under this dimensional reduction, the one-loop corrected metric of the AdS4 soliton geometry 
reads [33] 

dsl = f{r)de' + ^ + 72 (^4 + dx^), (30) 

g{r) L^ 



where 



5(r) = ^(l-j|).(l+7-/3(r)), 

1 / /f^o j>^o />->y \ 

a(r) = --^(l7^ + 17^-ll^), 
4L6 y r^ r6 r^ J 

/3(0 = -^ (17^ + 175- 67^ . (31) 



* Here we assume that the IID one-loop corrected action H25p may be dealt with classically on the spherical 
compactified background and on its dimensional reduction. We believe this is plausible because a sphere does not 
have noncontractable cycles, which would bring about new loop diagrams. The authors thank Shinji Hirano and 
Masaki Shigemori for discussions on this issue. 



The AdS radius L is related to the number N of M2-branes by 



and the volume of S" is given by 



^' = ^"'Wi^^' <=''' 



Vol(S') = y(2L)'. (33) 



The radius in the ^-direction is found to be 



2.i?=i^-^. (34) 

3ro SLVo ^ ' 



Using (j26p (j32p . the expansion parameter 7 is reexpressed as 



223/3.^16/3 ^6 

7= 5 1^- (35) 



It is useful to note 

Vol(5^ 



3 




N' 


1 




iV3/2 


23/2. 


37r 


L2 ■ 



(36) 
In the presence of the deficit angle 5 = 2tt{1 — n), the higher derivative term W4 behaves like 

^4 = 2 x4:i2 + 2vr(l - n)5(7A) • j^ + 0(7) + O ((1 - nf) , (37) 

using the uncorrected AdS4 soliton metric. This leads to the contribution 

Sr = ^g?'Vol(5'). (38) 

The contribution from the Einstein-Hilbert action can be found by applying the area formula to 
the corrected metric 

4-) . i^.v„„57r*f-5:i 



K 



11 Jro 



Att^R 7 ( ( '^1 



K 



•Vol(SO-(^oc-(l + ^)r-o) (39) 



Finally the total contribution to Sa is 

'^''^ ' ^4,^«,_^ + ?^ + 0(.)^|. (40) 



L2 23/2 . 3^ \^ °° 3 L4 

In the end, the 1/A^ correction to the finite term is found to be 

(-'^'" -'"").„.. - f ^^(-t::"'^""'^' 

= _^7V3/2 + !I^!^!^iVi/2 + o(Ar-i/2) 
9 9 -^ \ ) 



(41) 



2.2.2 M5-brane Soliton 

Finally we would like to study the AdSy soliton dual to the six dimensional SCFT of N M5-branes 
on a supersymmetry breaking circle. In a similar way to the AdS4 case, we reduce the IID theory 
on S'^ and then the one-loop corrected 7D metric reads [33] 

9v) ^^^ 

where 

/W = ^(l-;|J- (1 + 7 •«(-)), 

5W = ^(l-5l)-(l + 7-/3(r)), 

2-73 / r6 ri2 rl^\ 

The AdS radius L is related to the number N of M2-branes by 



and the volume of 5^ is given by 



L9 = iVslMi , (44) 



3 
The radius in the ^-direction is found to be 



Vol(5^) = ^{L/2f. (45) 



„ 27rL2 22.7-737r7 

2ttR = -, -. 46 

3ro 15L4ro ^ ' 



Using (I44p . the expansion parameter 7 (j26p is reexpressed as 



22/3 . ^16/3 ^6 

7 = 
It is useful to note 



^=^ Wr ^''> 



Yol{S^) _ 2 iV3 
" 3^" L5"' 



ni 



(48) 



In the presence of the deficit angle 5 = 27r(l — n), the higher derivative term VF7 behaves like 

Wr = ^s^r ° + Ml - n)6{^A) ■ ,y^,/° + 0(7) + O (d - n)^) , (49) 

using the uncorrected AdSy soliton metric. This leads to the contribution 

(VK) _ 2^ • 3 • 73^ r^ IttRVs • Vol(5^) 

^ ~ 5l6 14 2^^f; ^' 

73 2i4/3.7ri3/3 Ary., 

= 35 ^ V + OCV-')- (») 

9 



The contribution from the Einstein-Hilbert action can be found by applying the area formula to 
the corrected metric, as 

^^EH) ^ 47r^fiy3 • Vol(g^) r-^^r3 / 2^ • 13 • 73 • 7 r^\ 
SiV^yg^ /I ^ 1 ^ 2'^-13-73 7 



4 _ J^ 4 _ ^ xo 10 / 4 

SttLS V4''°° 4""° 5-14 'ls^" 

2iVV3i2 4 2 iVVs 2l''/3 . 73 . ^13/3 ]S[V3 

37rL8 ""^ ~ 357r ijs §^7^^ ^ 

Finally, the total sum of the finite term looks like 

2 N'-^V^ 

' finite S^TT R^ 

that is, the correction to Sa vanishes in the leading order (~ N^). 



+ 0{N-^). (51) 



K"' +<■•)_ - -i' 5^ -«(-'). (-) 



3 Holographic Entanglement Entropy in Gauss-Bonnet Gravities 

In this section, we will investigate holographic entanglement entropy for the half space of the 
boundary of AdS soliton geometries in Gauss-Bonnet gravities. 

3.1 4D Gauss-Bonnet Gravity 

We consider 4D AdS-Einstein-Gauss-Bonnet gravity|fj 



16^gS^^ 



-2K + n+-!-—CGB 



(53) 



where 

A = - -^ , Cgb = n^upan'^"'"' - ^n^un^'" + n^ . (54) 

The Euclidean AdS soliton metric (= Euclidean AdS-Schwarzchild metric) is 

ds^ = -^yJT^dr^ + "-^X^de^ + ^{dtl + dx^), (55) 

where 

/ ^ \ 3 2/2 

/(r) = l- - , 9^e + 2T:R, R=i^. (56) 

The solution does not depend on the Gauss-Bonnet coupling rj, since the Gauss-Bonnet term is 
topological in 4D. 



^ Precisely speaking, some surface terms should also be included in the action (except for a particular value of 
77 = 1/2 in the 4D case [34]) and they would affect the divergent part of the entanglement entropy. Since we focus 
on the finite part of it in this section, we simply omit them here. 

10 



We define the subsystem A of the boundary time shce t^; = by x > 0. To compute the 
entanglement entropy Sa holographicahy, we introduce a conical deficit on the (t_B,x)-plane along 
with the prescription of the replica method, and the bulk curvature tensors at that time is again 
given by the formula (J14p . Then the contributions of the Einstein-Hilbert term and Gauss-Bonnet 
term are computed as, respectively, 

d 



.{EH) 



y{EH) 



Area(7yi) 

4G 
2'kR 



n=l 



(4) 

N 



TtL^ Too 



dr 
7rL2 



3G 



(4) 

N 



and 



S 



(GB) 



^ 1 
on 



yiGB) 



Z- 



{GB)\n 
1 I 



2ttR 



32TrG[ 

TV 



27rL^ roc 



n=l 



dr^/g 



167r/(r) 



^•o 



3G 



2r 



^'o 



+ 



■ri + 



r / 

(4) 

N 



G' 



T] (roo -> oo). 



Therefore the total expression of Sa is, 

Sa = (div. 



7rL2 



3G 



+ 



^L2 



G 



V- 



3.2 5D Gauss-Bonnet Gravity 

Next we consider 5D AdS-Einstein-Gauss-Bonnet gravity. 



IGttG 



^/^'^^^ 



-2A + 7^ + 



7]L^' 



C-GB 



(57) 



(58) 



(59) 



(60) 



where A = — -^ and the expression of Cgb is the same as (|54p . The value of the Gauss-Bonnet 
coupling T] is restricted to the range of 



7 



9 



36 ' 100 
from the demand of causality in the Lorentzian version of this theory 



(61) 
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The Euclidean AdS soliton solution in this theory is given by the metric 



where 



ds' 



^« = ^^ 



L' 



r'^f{r 



-dr' + 






dr + 



^AdS 



{dt\ + dx^ + dy^), 



(62) 



'i-4^U-(7 



LacLS 



L 



To 



/oo = lim fir) 
0^6 + 2ttR, 


1 + VI - 4ry 



(63) 



It is useful to note that f'{ro) = 4/ro and f{ro) = linir-i^oo /'(^) = 0. This spacetime is asymptot- 
ically an AdS5 with the radius L^^s. This solution was found in 06] in the spherical form, as a 
black hole solutionis 

We define the subsystem ^4 as a; > again, and call the coordinate length along the y-direction 
Vy (we take Vy — )• oo). In this case. 



JEH) 



27rR ■ Vy '-''° 



AG 



(5) 

N 



ro ^AdS 



-dr 






Jo 



IQG%^R 



(64) 



and 



S 



{GB) 



2ttR-V, 



4G 



y V 

LAdS Jro 



dr [6/(r) + 6r/'(r) + r /"(r)] r 



^^|^^[3rV(r)+rY( 



ro 



2Gjy LAds 

STrVyRrl ,„ WyL^ r- 

'IJoo ' CK-i 'lyJo 



2GPl 



2G%^R 



The total expression of Sa is, 



^•(l-6./oo)/^ 

AGfl 



^"^^"^^ -(l-8.)v^ 



IQG%^R 



One may find that the divergent part looks proportional to the central charge 

''-"^^'^"-(l-ery/oo), 



8G 



(5) 

N 



(65) 



(66) 



(67) 



by an appropriate way of parameter fixinglj But we do not know whether it is sensible because 
the divergent part depends on the way of regularization. 



Before it, asymptotically flat black hole solution for A = had been found in [ST]- The planar form of the AdS 
one we are using here was presented in [SS] and discussed as an AdS soliton in [39j . 

"^ Fixing LacIS, R, Vy, and Roo = Rtoo/LacIS realizes it, where Roo is the "physical radius" of the 6'-circle on the 
boundary. 
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4 Entropic Phase Transitions in 5D Gauss-Bonnet Gravity 

Here we would like to calculate the holographic entanglement entropy (HEE) for a subsystem A 
which is defined by a strip with the width / in Gauss-Bonnet theories. We are interested in a 
discontinuous transition on the behavior of the holographic entanglement entropy for the AdS 
soliton as we change the size /. In the absence of the higher derivatives, this has been studied 
in |23H25] and this phase transition is identified with the confinement/deconfinement transition. 
There is a critical value I = Ic such that -^ = (confined phase) for / > Ic and -^ > (deconfined 
phase) for / < Ic- Qualitatively similar results has been recently reproduced in lattice gauge theory 
approaches [26H28] . Below we would like to study how this behavior is affected by higher derivative 
terms. 

In the four dimensional Gauss-Bonnet gravity, the curvature correction term in the HEE formula 
(l2|) is topological and is determined only the Euler number of the surface 7a • Thus the I dependence 
of its higher derivative correction gets trivial. Therefore below we will study the 5D Gauss-Bonnet 
gravity. 

4.1 Extremum equation for holographic entanglement entropy 

We employ the proposed generalization ([2]) |161ll7j of the original HEE [2] in the Gauss-Bonnet 
gravity. More precisely, including the boundary term, the HEE is given by minimizing the functional 

Sa = — J- dx^Vh{l + r]L'^n) + ^^^^ dx^^/hlC, (68) 

where h is the induced metric on the three dimensional surface 'ja which satisfies d^A = dA and 
TZ is the intrinsic curvature of h; hi^ is the induced metric on djA and /C is the trace of its extrinsic 
curvature. If we define a unit normal vector of the embedding of d'jA in 'J A by n", then we have 

dx'^^/hb-^= I dx^da{^/h'b-n'')\r=r^. (69) 

We assume that the subsystem A sits on the time slice i = and extends in y- and ^-directions. 
Therefore we can specify the profile of ja by the embedding function 

X = x{r). (70) 

Its holographic dual surface 7a is given by a codimension three surface defined by 

t = 0, x = x{r), -oo<y <oo, 0<e < 2-nR. (71) 

The regularized length in y is denoted by Vy below. The surface 7a is determined by extremizing 
the functional (|68p . When there are several extremal surfaces, we select the one which has the 

13 



smallest value of (|68|) . We choose the boundary condition of 7^ such that it starts from x = ^ 
at r = Too, extends into the smaller r region until it reaches r = r^^, and comes back to the AdS 
boundary r = Too at a; = — 2- 

The induced metric on 7^ is given by 

ds^ = Kbdx-dx' = (^^ + -^{x'{r))A dr' + "^^^dO^ + -^dy\ (72) 



We obtain 



L-L2 









,M = _^^m±£m=. (74) 



where we defined 

q{r) = - 

Ly/L' ■ L\,, + r^f{r)x>^ 

When we integrate on 7^1, the second term q'{r) in (|73p yields surface terms. We can show that the 
contribution of q{roo) is canceled by the Gibbons-Hawking term in (|68|) . Furthermore, since the 
profile is connected with its mirror image at r = r,, , there is indeed no boundary there and neither 
is the boundary term g(r^). Therefore the functional we need to minimize is given by 

S^ = ^^^ r^ L^-L\^s-^i^ + 9{r))+r^f{r)x'^ ^^ 



where 



TiRVy f- rg{r) + rX_^^ 



G^^ ■ LAdS Jr, ^ 



g{r) = 2r]f{r)+r]rf'{r), (76) 

^ = 1 + 72^rV(r)x'^ (77) 



and r = r* is the turning point of the surface. Note that X > 1. 
The equation of motion for (|75p leads to the conservation law 

r>f{r)x'{L^ ■ L\,, + r'f{r)x'' - L^L\,,g{r)) _ 3 



VTR, (78) 



where r^ is some constant. If the surface is smooth at r = r*, i.e., x'(r*) gets divergent, then 
rs = r* . We can find the width of the subsystem A from 

rroo 

1 = 2 dr x'{r). (79) 
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Now let us examine the equation ([75]) . We impose the obvious inequahty: 

/oo>0 ^ v<\- (80) 

Although the causality bounds for i] is known to be given by (j6ip [35], at first we would like to 
proceed by ignoring it and would like to come back to this in the end. By taking the square of 
([78|) . we obtain 

where 



0, (81) 



F{X) = (rV(r)-rf/(r,))X3-rV(r)(l + 2<7(r))X2 

+r6/(r)(25(r) + g{rf)X - r^f{r)g{rf. (82) 

4.2 Various phases for "minimal surfaces" 

Let us examine the solutions of the equation (j78p . In order to have a physical solution for the 
surface 7^, we need a solution to the cubic equation (j8ip with X{r) > 1. This constrains the range 
of r where desirable solutions exist. 

The existence of smooth solutions to (j8ip requires that X is positively divergent at the turning 
point r = r^. In this case we find r* = r^ as follows from (j78p . This requires that a solution to 
()8ip gets positively divergent as we approach r = rg from the above. This condition leads to the 
inequality g{rs) > —1/2 or equally, 



1 ff - 9 + 3v^9ff - 2f 4 + 9 
32 Vi 



V>-:^- ^ ' , (83) 



s 



where f^ = rg/ro. This is always satisfied when 



Also, for the values 



only a particular range of r^ 



V > --• (84) 



To < " * -8' '-''^ 



is allowed. For ^ < — ^ there will be no smooth solutions. See Fig (1] for numerical plots of the 
bound. 
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rl = 9/100 

II = -1/8 

11 = -7/36 

11 = -5/16 




Figure 1: The bounds for g{rs) > —1/2 (|83p(|86p (red) on the (rs/ro,?7)-plane. Lines for rj = 
1/4, 9/100, -1/8, -7/36, -5/16 are also displayed. 

In these cases, a solution to (j8ip satisfies X > 1 for r > rg and it diverges at r = r^ while there 
is no such solution for r < Vg- Moreover by studying the second order perturbations, it is indeed 
a local minimum of the functional (j68p . In this way, we find a smooth solution with the correct 
boundary condition and thus this gives the first candidate for ^a (as depicted in Figl2]|^a)). 

Now we would like to turn to the other cases with g{rs) < —1/2. Since we cannot find smooth 
solutions, we need to be satisfied with a solution which is not smooth at the turning point r = r^, as 
e.g. the one depicted in Figl2]^b), assuming that the minimal surface always existqj. The existence 
of solution leads to another bound r > r^, where r^ satisfies 



r!f{rs 



4 (gin) - 1)3 



■rifirb). 



(87) 



27 g{n) 

For r > Tf,, there are two solutions which satisfy X > 1. The smaller is continuous through 
r = rs{> Tfi), and the greater is divergent in the limit r — t- r^ — 0. They get coincident at r = r^, 
and disappear for r < rf,. Joining these, we can consider, for example, a "switch-back" surface like 
FigEI^c). However, the second order perturbation of (|75p is positive for the smaller solution, and 
negative for the greater. Then we need to consider only the first, and in this case the form of the 
surface is always like Figj2]^b). 

Other than the solutions discussed above, there are two "trivial" candidates of the "minimal 
surface". One is displayed in Figl3|^d), which is a pair of the surfaces x = —1/2 and x = 1/2 
connected at r = tq. It is not a solution of the extremum equation F[X) = but on the edge of 



*In the previous case with g(rs) > —1/2, we can also construct a similar non-smooth surface like Fig[2jb). However, 
we can confirm that the functional (|68p is always greater than the smooth one. 
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(a) smooth surface (b) surface singular at r = r* (c) "switch-back" surface 

Figure 2: Patterns of extremal curved surfaces ^a (projected on the (r, a;)-plane). The integrand 
of ()75p is locally minimum on the red lines and maximum on the blue one. Therefore we only have 
to consider (a) and (b). 



(d) connected at r = tq (e) disconnected 

Figure 3: Two patterns of "trivial" ^a- They have different values of Sa due to the boundary 
term g(ro), accounted only for (e). For r/ < 0, (e) is preferred to (d), and vice versa for r/ > 0. 
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the configuration space of x{r). The corresponding value of Sa is, independently from /, given by 
substituting r* = vq and X = 1 to ([75]) . resulting 



Sf = l^^y r{rg{r)+r)dr 



2GjY • LAds 

The other, which is displayed in Figl3]^e), is similar to it but different in an important way. Unlike 
the all variations of the surfaces discussed above, it consists of two parts which are not connected to 
each other. Therefore the corresponding value of Sa includes the contribution from the boundary 
term g(ro) = -Arl/L'^LAdS and 

q{e) _ ^(d) 2TTRVyr]L'^ 

^A - ^A —(5) lyO) 

- "^^' [il + 2rjUrl-il-8v)rl]. (89) 



2Gj^ ■ LAds 



This is, as is expected, precisely twice of (166p except the difference of the divergent part coming 
from the boundary termo We can immediately see that, when ry < the phase (e) is preferred to 
(d), and for r/ > (d) is preferred. 

4.3 Numerical results and phase structures 

Now we can numerically execute the integrals (I75p and ()79p for the solutions of the extremum 
equation (I78p . and plot the values on the {I, /S'yi)-plane. We set the turning point r* as r* = r^ when 
9{i~s) > — ^, while we set it as r* = ri, when g{rs) < — ^. These are because we find that they give 
the smallest values of the functional (j68p when we change r^, with I fixed. The resulting behavior 
of Sa{^) is plotted in FigH] (?/ > 0) and Figj5] (?/ < 0). The zero point of Sa is taken to be S)^ 
()88p for ?7 > and S^ ()89p for 77 < 0. The case for r/ = is the same one that was investigated 
in [231I21]- 

Where the curve Sa = Sa{1) is below the /-axis, the corresponding nontrivially connected 
surface, given as the solution of the extremum equation, is the "minimal surface" 7^ for each I. 
Otherwise the minimal surface is given by the "trivial" surfaces (d) or (e), depending on the sign 
of r]. Therefore, when the line crosses the /-axis, there occurs a phase transition between those 
different phases as we mentioned before. 



" In this section we take the last term in (|68p into consideration, while we omitted the boundary term in (|6U|I . In 
fact, if we were to omit the boundary term in (|68p . we have an additional contribution from q{r 00) to (|89|l and the 



result agrees with twice of ([66} exactly. 
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Figure 4: {1,Sa) plots for variable < i] < 1/4. The unit for I and Sa are R and TTVyL/2G}^ , 
respectively. The value for the phase (d) is taken to be for Sa- 



First look at FigUJ We can see that the qualitative form of the curve does not change in the 
range < 77 < 1/4. The phase of the curved surface is preferred when / is small, and there occurs 
a phase transition at some particular value / = Ic, for example, Ic/R — 0.338 for 1] = 0.05. Notice 
that nothing special happens around the causality upper bound 77 = 0.09 (I6ip . 

The plots for r/ < 0, displayed in FigUl may be more interesting. As we lower the value of 77, 
the plotting curve changes its shape. Below r] = —1/8 = —0.125, the point corresponding to the 
boundary r* = Tc (j86p between (a) and (b) walks on the curve from the origin toward the turning 
point. Just before it reaches there, when rj ~ —0.209, one (or more) loop(s) appears around there 
(FiglH]). After that Tc continues to walk downward until i] = —5/16 = —0.3125. Regardless the 
phenomena above, the phase structure is not altered. However, across rj = —0.548 . . . , the shape of 
the plotting curve changes dramatically, and below there, the curve is above the l-axis everywhere 
and so the disconnected phase is always favored (Fig. [7j). Note that these strange phenomena do 
not take places above the causality lower bound rj = —7/36 = —0.1944. . . (I6ip . 

In summary, if we impose the causality bound (I6ip . none of new phenomena found for the higher 
derivative theory such as the non-smooth surface as ja like Figl2|^b), the loopy profile of Sa{1), 
or the absence of phase transition, do not occur. In this sense, we can conclude that the higher 
derivative corrections does not largely change qualitative properties of HEE for AdS solitons. If we 
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Figure 5: {1,Sa) plots for variable r^ < 0. The value for the disconnected phase (e) is taken to be 
for Sa- 
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temporally forget the causality bound, then from the above analysis we can learn that there should 
be lower bound for rj otherwise the strange behaviors start to occur. For example, if we require 
that the non-smooth surfaces should not appear as 7a, then we obtain the bound 7] > —5/16. Also, 
if we exclude the absence of phase transition, then we find rj > —0.548 • ••. Moreover, if we consider 
the Gauss-Bonnet gravity literally without considering more higher derivative terms, then we will 
find that HEE is not well-defined for the region ry > as we will explain in the next subsection. 



4.4 Comments on an instability for i] > 

Before we conclude this paper, we would like to point out an important fact for the conjectural 
formula ()68p in the Gauss-Bonnet gravity. Consider a smooth surface 7^ which minimizes the 
functional (j68p . We can assume that it is symmetric along 6. Let us focus on the region near the 
turning point r = r^, where we can treat the warp factor of the AdS space as a constant. In that 
region, we can simply ignore the ^-direction and we can regard ^a as an effectively 2-dimensional 
surface in a flat 3-dimensional ambient space, spanned by {r,x,y). In this setup, we can add 
infinitesimally small handles to ja near r = r^: and decrease the value of the f\/hTZ (i.e. the Euler 
number) without changing the other terms. Thus if we assume r/ > 0, we can take (j68p to be 
an arbitrarily small value and so there is no minimum of the this functional. In this way we find 
that HEE Sa is ill-defined for 77 < in the Gauss-Bonnet gravity. This argument gets more clear 
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in 4D Gauss-Bonnet gravity as the curvature contribution in ([55|) gets purely topological without 
focusing on near the turning points. 

One may worry that our analysis for r/ > may be meaningless. However, this is not the case 
if we implicitly assume the presence of more higher derivative corrections in addition to the Gauss- 
Bonnet gravity, which will be the case in string theory. Indeed, we can show that the problem we 
mentioned occurs only when the higher order corrections are absent. In other words, what we find 
here is that if we consider the purely Gauss-Bonnet gravity in any dimension literally without any 
more higher derivative terms, then we will find that HEE is not well-defined when r/ > 0. 

5 Discussions and Summary 

In this paper, we studied the holographic entanglement entropy (HEE) for AdS soliton geometries 
in the presence of higher derivative corrections. Our results in this paper show that the proposed 
higher derivative correction to HEE (|68p behaves in a sensible way. 

In the first half part, we calculated the leading higher derivative corrections due to TZ'^ term 
for AdS soliton geometries in string theory and M-theory. Via the AdS/CFT, they are dual to 
the strong coupling expansions in the dual confining gauge theories such as a compactified AA = 4 
super Yang-Mills. Our result is qualitatively consistent with the free Yang-Mills calculation of the 
entanglement entropy. 

In the latter half part, we studied the HEE for AdS solitons in the Gauss-Bonnet gravity. 
Especially we examined the dependence of HEE on the size of the subsystem A. If we restrict to 
the known causality bound of the Gauss-Bonnet parameter ?/, our result shows that the behavior 
of HEE and the structure of phase transition are qualitatively similar to that for Einstein gravity 
i.e. 7/ = 0. However, if we go beyond that bound, we observed several strange behavior such as the 
absence of phase transition and singular behavior of the solutions. We also find that if we consider 
the purely Gauss-Bonnet gravity literally without any more higher derivative terms, then we will 
find that HEE is not well-defined when r/ > 0. 
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